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ABSTRACT 

Aims. This research has two main goals. First, we present the atmospheric structure and the fundamental parameters of three red 
supergiants (RSGs), increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible 
mechanisms that may explain the large observed atmospheric extensions of RSGs. 

Methods. We carried out spectro-interferometric observations of the RSGs V602 Car, HD 95687, and HD 183589 in the near-infrared 
A'-band (1.92-2.47 pm) with the VLTI/AMBER instrument at medium spectral resolution (R —15001. To categorize and comprehend 
the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D 
convection, and new 1-D self-excited pulsation models of RSGs. 

Results. Our near-infrared flux spectra of V602 Car, HD 95687, and HD 183589 are well reproduced by the PHOENIX model atmo¬ 
spheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to 
determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, 
the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. 
Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric 
structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed 
decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing 
observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temper¬ 
ature or variability amplitude. 

Conclusions. The location of our RSG sources in the Hertzsprung-Russell diagram is confirmed to be consistent with the red limits of 
recent evolutionary tracks. The observed extensions of the atmospheric layers of our sample of RSGs are comparable to those of Mira 
stars. This phenomenon is not predicted by any of the considered model atmospheres including available 3-D convection and new 1-D 
pulsation models of RSGs. This confirms that neither convection nor pulsation alone can levitate the molecular atmospheres of RSGs. 
Our observed correlation of atmospheric extension with luminosity supports a scenario of radiative acceleration on Doppler-shifted 
molecular lines. 

Key words, supergiants - Star: fundamental parameters - Star: atmospheres - Star: individual: V602 Car. HD 95687 and HD 183589. 


1. Introduction 

Red supergiant (RSG) stars are known to lose mass with mass- 
loss rates of 2 x 10 7 M e /yr - 3 x 10 4 M Q /yr (de Beck et al. 
12010 ). and they are one of the major sources of the chemical 
enrichment of galaxies and of dust in the universe, along with 
asymptotic giant branch (AGB) stars and supernovae. Currently, 
the mechanism of mass loss of RSG stars and semi-regular or ir¬ 
regular AGBs is not known in detail. Nevertheless, in the case 


* Based on observations made with the VLT Interferometer (VLTI) 
at Paranal Observatory under programme ID 091.D-0275 


of Mira-variable AGB stars (mass-loss rates of 10 6 M Q /yr - 
1CU 4 M 0 /yr, Wood et al. 119831 19921 , this mechanism is better 
understood. The theoretical models that explain the Mira mass- 
loss process are based on pulsations that extend the atmospheres 
to radii where dust can form, and subsequently on radiative pres¬ 
sure on dust grains that drives the wind (e.g., Bladh et al. 120131 . 
In the case of variable RSGs, the amplitude of the light curves is 
about one-third of that of Miras (e.g.. Wood et al. 119831 ), so that 
pulsation is expected to play a less dominant role (cf. Josselin 
& Plez 12007 ). Other mechanisms that might give rise to mass 
loss in the RSGs are convection and rotation (e.g., Langer & 
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Heger !1999 ). Constraints on the mechanisms that levitate the at¬ 
mospheres of RSGs are thus fundamental for our understanding 
of the mass-loss process of RSG stars. 

Observations of Mira variable stars using the IOTA or VLTI 
interferometers show evidence of molecular layers lying above 
the photospheric layers (e.g., Perrin et al. 1 2004 ; Wittkowski et al. 
120081 ' 2011 1. Theoretical dynamic model atmospheres (Ireland 
et al. 12004a, 2 004bl 120081120TT1 Scholz et al. 120131 can explain 
reasonably well these molecular layers for Miras. On the other 
hand, interferometric observations of RSGs also indicate the 
presence of extended molecular layers (CO and water), which 
cannot be explained by hydrostatic model atmospheres (Perrin 
et al. 120051 Ohnaka et al. EoTTl I20l3l Wittko wski et al. I20l2l >. 
The red supergiant VX Sgr (Chiavassa et al. 12010al ) showed a 
good agreement with Mira models, although they have very dif¬ 
ferent stellar parameters than expected for this source. 

This paper is conceived as part of a series of three previous 
papers (Wittkowski et al. l20l2 r ; Arroyo-Torres et al. l2013lf2014l 
In these previous works, we presented the atmosphere structure 
and the fundamental parameters of a sample of four RSG stars 
and five cool giant stars. One goal of the current paper is to add 
three red supergiants and thus to increase this sample. On the 
other hand, our previous works showed that the observed visibil¬ 
ity data of the RSGs and of one of the red giants, /3 Peg, indicate 
large extensions of the molecular layers, similar as those previ¬ 
ously observed for Mira variable stars (Wittkowski et al. [20081 
12011 ). This was not predicted by hydrostatic PHOENIX model at¬ 
mospheres. However, the spectra of all our stars were reproduced 
well by the PHOENIX models. This indicates that the molecular 
opacities were adequately included in these model atmospheres, 
but that they were too compact compared to observations. In or¬ 
der to understand the processes that may explain the extended 
molecular layers, the second goal of this paper is to investigate 
the effects of realistic three-dimensional (3-D) radiative hydro- 
dynamical (RHD) simulations of stellar convection as well as of 
one-dimensional (1-D) self-excited pulsation models on the ex¬ 
tensions of RSG atmospheres. These processes were previously 
discussed as possible mechanisms to levitate RSG atmospheres 
(Chiavassa et al. l2010a|l2011a i). 

In this paper, in addition to studying RSGs, we also refer to 
asymptotic giant branch stars (AGBs), Mira variable stars, and 
red giant stars. AGB stars are low and intermediate mass evolved 
stars before they evolve toward hotter temperatures in the HR di¬ 
agram. Mira stars are long-period large-amplitude variable AGB 
stars. With red giants we refer to giants on the first red giant 
branch. 

Our work is structured as follows: In Sect. [2] we describe 
our AMBER observations and the data reduction. In Sect. [3] we 
present the results obtained from the PHOENIX model fitting 
and the fundamental parameters. In Sect. [4] we characterize the 
extensions of the atmospheres. In Sect. 0 we show the results 
obtained from the comparison with the convection and pulsa¬ 
tion models and we discuss alternative mechanisms. Finally, in 
Sect. [6] we summarize our results and conclusions. 


2. Observations and data reduction 

We observed V602 Car (Simbad spectral type M3-M4 I), 
HD 95687 (M3 lab), and HD 183589 (K5 lb) with the ESO 
Very Large Telescope Interferometer (VLTI), utilizing three of 
the Auxiliary Telescopes of 1.8 m diameter, and the AMBER 
instrument (Astronomical Multi-BEam combineR) with the ex¬ 
ternal fringe tracker FINITO (Petrov et al. [2007). We worked 


in medium-resolution mode (R ~ 1500) in the K-2.1 pm and K- 
2.3 pm bands. The integration time (DIT) of each frame was 20 
ms. Our data were observed as sequences of cal-sci-cal (cal is 
calibrator and sci is our target), with 5 scans for each of them. 
Table [T|lists the detailed information about our observations and 
the calibrator used for each target. Table[2]shows the calibrators 
used for our observations together with their angular diameters. 
We selected them from the ESO Calibration Selector CalVin, in 
turn based on the catalog of Lafrasse et al. (120101 ). 

During the acquisition of one AMBER frame, there are op¬ 
tical path fluctuations (jitter) that produce fringe motions. These 
motions reduce the squared visibility by a factor e FINITO 
factor, where <x (t is the fringe phase standard deviation over 
the frame acquisition time. This attenuation is corrected in the 
science data by the calibration provided the FINITO factors 
are similar in science and calibrator (more information in the 
AMBER User Manual []). 

As a first step we selected the scans such that the FINITO 
factors were similar between cal and sci data. After that, we ob¬ 
tained the visibility data from our selected AMBER observations 
using the 3.0.7 version of the amdlib data reduction package 
(Tatulli et al. 120071 Chelli et al. 120091 ) . This included the removal 
of the bad pixel map and the correction for the flat contribution. 
Afterwards, we calculated the pixel-to-visibility matrix (P2VM) 
to calibrate our data for the instrumental dispersive effects, and 
obtained the interferometric observables (visibility and closure 
phase). Next, we appended all scans of the same source taken 
consecutively, selected and averaged the resulting visibilities of 
each frame using appropriate criteria. In our case, the criteria 
were based on the flux (we selected all frames having flux densi¬ 
ties three times higher than the noise) and on the signal-to-noise 
ratio (S/N). We only used 80% of the remaining frames with best 
S/N0. 

Using scripts of IDL (Interactive Data Language), which 
have been developed by us, we performed the absolute wave¬ 
length calibration by correlating the AMBER flux spectra with 
a reference spectrum, that of the star BS 4432 (spectral type 
K4.5 III, similar to our calibrators; Lan£on & Wood 2000 ). A 
relative flux calibration of the targets was performed by using 
the instrumental response, estimated by the calibrators and the 
BS 4432 spectrum. Finally, calibrated visibility spectra were ob¬ 
tained by using the average of two transfer function measure¬ 
ments taken before and after each science target observation. In 
the case of V602 Car (25-26 April), HD 95687 (24-25 April) 
and HD 183589 (03-04 August), we only used the first calibra¬ 
tor, because the other calibrator had very different FINITO fac¬ 
tors and their visibilities were not of sufficient quality. The error 
of the transfer function was calculated as in our previous work 
(Arroyo-Torres et al. 1201 31 120141) . 

3. Results 

We compared our observational data to synthetic data pro¬ 
vided by a grid of PHOENIX model atmospheres (version 16.03, 
Hauschildt & Baron [19991 from Arroyo-Torres et al. 20131) . 
These models are based on a hydrostatic atmosphere, local 
thermodynamic equilibrium and spherical geometry. The best- 
fit PHOENIX model is obtained from an iterative process using 
the continuum band around 2.25 pm, as explained by Arroyo- 
Torres et al. (12013) . The final values for the used PHOENIX 

1 http://www.eso.org/sci/facilities/paranal/instruments/amber/doc.html 

2 see AMBER Data Reduction Software User Manual; 
http://www.jmmc.fr/doc/approved/JMMC-MAN-2720-0001 .pdf 
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Table 1. VLTI/AMBER observations 


Target (Sp. type) 

Date 

2013- 

Mode 

K- (urn) 

Baseline 

Projected baseline 
m 

PA 

deg 

Calibrator 


04-04 

2.1 

A1-G1-K0 

75.86/80.02/127.8 

91/21/55 


V602 Car (M3-M4 I) 

04-04 

2.3 

A1-G1-K0 

79.89/68.52/112.6 

130/48/93 

HR 4164-z Car 


04-26 

2.3 

D0-H0-G1 

63.22/58.73/64.46 

71/-172/125 



04-04 

2.1 

A1-G1-K0 

73.06/80.57/128.6 

80/14/45 


HD95687 (M3 lab) 

04-04 

2.3 

A1-G1-K0 

79.94/71.35/117.5 

121/43/85 

HR 4164-z Car 


04-25 

2.3 

D0-H0-G1 

56.88/53.89/69.30 

104/-153/153 



05-04 

2.1 

D0-I1-G1 

78.89/45.56/63.26 

99/-134/134 


HD 183589 (K5 lb) 

07-29 

2.3 

D0-I1-G1 

69.75/46.61/55.98 

99/-134/141 

38 Aql - HR 7404 


08-04 

2.1 

A1-G1-K0 

88.87/67.14/126.2 

-145/-72/-114 



Notes. Details of our observations. The projected baseline is the projected baseline length for the AT VLTI baseline used, and PA is the position 
angle of the baseline (north through east). 
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Fig. 1. Normalized flux (a), squared visibility amplitudes (b), UD diameters (c), and closure phases (d) for the example of V602 Car 
obtained with the MR-K 2.3 /u m setting, on 26 April 2013. In black the observed data, in blue the UD model, and in red the PHOENIX 
model. 


model atmosphere are: For V602 Car, T e ff=3400 K, log(g)=-0.5; 
for HD 95687, T eff =3400K, log(g)=0.0; and for HD 183589, 
T e ff=3700K, log(g)=1.0. For all cases, we used a model with 
solar metallicity and a micro-turbulent velocity of 2km/s. We 
chose a mass of 20 M e for V602 Car and HD 95687 and of 1 M e 
for HD 183589. We chose a low mass of 1 M Q for the latter tar¬ 
get, because the final parameters indicate that it is a source with 
lower luminosity and thus lower mass compared to the other 
RSG sources. We note that the structure of the atmosphere is not 
very sensitive to variations of the mass (Hauschildt et al. 119991) . 
Certainly, any of those structure variations are below the level of 
the detectability of our interferometer. 

Fig. |U shows as an example the resulting normalized flux, 
squared visibility amplitude, uniform disk diameter, and closure 
phase data for one of our sources, V602 Car, obtained on 26 
April 2013. Also shown are the best-fit uniform disk model (blue 


curve) and the best-fit PHOENIX model atmosphere (red curve). 
The data and best-fit models for the remaining data are shown in 
the online appendix (Figs. [2^[6]). 

The normalized flux spectra show typical spectra of red su¬ 
pergiants in the K band (cf. Langon et al. 12007 ; Arroyo-Torres et 
al. 12013 ). The flux variations at wavelengths below about 2.0/rni 
are due to a higher noise level, possibly caused by the lower 
atmospheric transmission. In the K-2.3 band, we observe a de¬ 
creasing flux longwards of 2.25 fim and strong absorption lines 
of CO. The synthetic spectra of the PHOENIX model atmosphere 
are in a good agreement with our flux spectra including the CO 
bandheads. This indicates that the opacities of CO are well re¬ 
produced by the PHOENIX model atmosphere. 

The continuum visibility values near 2.25 n m are consistent 
with the predictions by the PHOENIX model atmospheres for all 
our sources. In the case of HD 183589, the visibility spectrum is 
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Table 2. Calibration sources 



Spectral type 

Angular diameter (mas) 

HR 4164 

K1 III 

1.64±0.12 

z Car 

M6 

1.54±0.11 

38 Aql 

K3 III 

2.22±0.02 

HR 7404 

K2 

1.17±0.08 


featureless and consistent with the PHOENIX model atmosphere 
prediction across the whole observed wavelength range. In par¬ 
ticular, the visibility spectrum of this source does not show fea¬ 
tures at the locations of the CO bandheads, which are visible in 
the flux spectrum, indicating a compact atmospheric structure 
where the CO layers are located close to the continuum-forming 
layers. Nonetheless, V602 Car and HD 95687 show large drops 
of the visibility in the CO bandheads between 2.3 /tm and 2.5 yum 
that are not reproduced by the PHOENIX model atmosphere. The 
synthetic PHOENIX visibility spectra show features in the CO 
lines, but these are much weaker than the observed features. 
This effect is also reflected in the panels showing the uniform 
disk diameter. The size increases of UD fits at the CO bandheads 
are about 40% for V602 Car and 20% for HD 95687, while the 
PHOENIX models predict UD size increases below 5%. These re¬ 
sults indicate that these sources exhibit a large contribution from 
extended atmospheric layers in the CO bands. The PHOENIX 
model structures are too compact compared to our observations 
for these two sources. We also observe a monotonic decrease be¬ 
yond 2.3 yum, which may be caused by pseudo-continuum con¬ 
tributions from CO or by contributions from water vapor. We 
observed the same phenomenon previously for the red super¬ 
giants VY CMa (Wittkowski et al. [20T2t i. AH Sco, UY Set, and 
KW Sgr (Arroyo-Torres et al. (120131 1. as well as for the small- 
amplitude pulsating red giants RS Cap (Marti-Vidal et al. 120 1 1 b . 
BK Vir (Ohnaka et al. l2UT2t . a Tau (Ohnaka l20l3bl ). and [3 Peg 
(Arroyo-Torres et al. 120141 ). 

The closure phase data of our sources in the 4th panels of 
Fig. |T|and Figs. [2]-[6]show variations within the noise level, and 
are thus are not indicative of deviations from point symmetry. 
However, since our measurements lie in the first visibility lobe 
and the noise level is relatively high, we cannot exclude asym¬ 
metries on scales smaller than the observed stellar disk. We note 
that there are points in the observed closure phases whose de¬ 
viation from zero is larger than the error bars. In general, small 
deviations from zero closure phases might indicate asymmetries 
in layers corresponding to certain atomic or molecular bands 
as previously observed for RSGs by e.g., Ohnaka et al. (1201 Ik 
Wittkowski et al. (120121 ). Ohnaka et al. (120131 ). However, it is not 
clear whether in our case this deviation are real or whether they 
correspond to systematic uncertainties of the data reduction, as 
for instance due to the bad pixel mask. 

3.1. Estimate of the angular diameter 

The continuum band near 2.25gm appears to be largely free of 
contaminations by molecular layers. Thus, fits of PHOENIX mod¬ 
els to the continuum band allow us to estimate reliable angular 
diameters of our sources. The angular diameter, obtained in this 
way, corresponds to the size of the outermost model layer (0% 
intensity radius). To estimate the Rosseland angular diameter 
(corresponding to the layer where the Rosseland optical depth 
equals 2/3), we multiplied our value of the angular diameter by 
the ratio between the Rosseland layer and the outermost model 
layer. This ratio was 0.92 for V602 Car, 0.95 for HD 95687, 
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Fig. 7. Squared visibility amplitudes in the continuum bandpass 
for V602 Car, HD 95687, and HD 183589 (from top to bot¬ 
tom) as a function of spatial frequency. Each point represents an 
average of data points within the continuum bandpass at 2.15- 
2.25 /rm. Shown are data of all observing dates and both spectral 
setups. The red lines indicate the best-fit UD models and the blue 
lines (often indistinguishable from the red) the best-fit PHOENIX 
models. The dashed lines indicate the maximum and minimum 
visibility curves, from which we estimated the angular diameter 
errors. 


and 0.93 for HD 183589. The model fits used all available data 
taken during all nights and with any of the two spectral setups, 
as both setups include the continuum band near 2.25 yum. Tab. [3 
lists the resulting best-fit Rosseland angular diameters as well as 
the best-fit UD diameters. Fig. 0 shows the continuum visibility 
data as a function of spatial frequency together with the best-fit 
PHOENIX and UD models. The errors of the continuum visibili¬ 
ties data were computed as an average of the individual errors, 
whereas, the errors of the angular diameter are estimated from 
the differences between the visibility curves lying at the maxi¬ 
mum and minimum of our data as shown by the dashed lines in 
Rg. m Deviating visibility points are caused by remaining sys¬ 
tematic uncertainties of the absolute visibility calibration. The 


4 































B. Arroyo-Torres et al.: Large extensions of RSG atmospheres 


Table 3. Summary of estimated angular diameters 



V602 Car 

HD 95687 

HD 183589 

0ud (mas) 
0Ro SS (mas) 

4.94±0.75 

5.08±0.75 

3.17±0.50 

3.26±0.50 

2.95±0.50 

3.04±0.50 
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Fig. 8. Effective temperature versus spectral type of our sources, 
together with calibrations of the effective temperature scale by 
Dyck et al. d 19981 ), Levesque et al. ( 120051 . and van Belle et al. 
( 120091 . Also included are previous measurements of RSGs and 
red giants as listed in the main text. In blue are the RSGs and in 
red the red giants. 


data of V602 Car include two points near the first visibility null, 
which increases the precision of the best-fit angular diameter. 

3.2. Fundamental parameters 

We estimated the fundamental stellar parameters of our sources 
to place them on the HR diagram. In particular we calculated 
the effective temperature, the luminosity, the Rosseland-mean 
radius, the bolometric flux, and the distance in the same way 
as described in detail by Arroyo-Torres et al. (120141) . We used 
the BVJHKs magnitudes from Kharchenko (120011 and Cutri et 
al. (120031 and the IRAS flux from IRAS ( 119881 . To convert the 
magnitudes into fluxes, we used the zero values from Johnson 
d 19651 and Cohen et al. (120031 . To deredden the flux values we 
used the color excess method applied to (V-K) and based on in¬ 
trinsic colors from Ducati et al. d20Qll . as described in Arroyo- 
Torres et al. d20141 . 

V602 Car and HD 95687 belong to the cluster CAR OB2 and 
we use the distance as determined by Humphreys et al. (119781 , 
For HD 183589, we used the distance value from van Leeuwen 
d20071 . Lastly, the effective temperature is estimated from the 
angular diameter and the bolometric flux, the luminosity from 
the bolometric flux and the distance, and the Rosseland radius 
from the Rosseland angular diameter and the distance. We as¬ 
sumed a 15% error in the flux, and a 10% error in the distance 
for HD 183589. For V602 Car and HD 95687, we used the er¬ 
rors from Humphreys et al. d!978b . The errors in the luminosity, 
effective temperature and radius were estimated by error prop¬ 
agation. The resulting fundamental parameters and their errors 
are listed in Table [4] 

The radius and luminosity of HD 183589 suggest that this is 
a source of lower luminosity and thus lower mass compared to 
the other observed RSG sources. It is at the limit between RSG 
and super-AGB stars. The visibility functions resemble those of 
the red giants as observed by Arroyo-Torres et al. (120141 . which 
do not show indications of an extended molecular layer. 



Fig. 9. Location of our sources in the HR diagram, compared 
to evolutionary tracks from Ekstrom et al. (120121) for masses of 
3 Mq, 4 M e , 5 M q , 7 M q , 9 M 0 , 12 M 0 , 15 M 0 , 20 M 0 , 25 M 0 , 
32 Mq, and 40 M 0 . The solid lines indicate models without ro¬ 
tation, and the dashed lines with rotation. Also shown are previ¬ 
ously measured sources as listed in the main text. In blue are the 
RSG stars, and in red the red giants. 


In Fig. [8] we plot the resulting effective temperatures vs. 
spectral types of our targets. For comparison, we include the 
calibrations of the effective temperature scale by Dyck et al. 
(1 19981 ) for cool giants stars, by van Belle et al. (120091 ) for cool 
giants stars and RSG stars, and by Levesque et al. (12005b for 
only RSGs. Fig. [9] shows the position of our targets in the HR 
diagram, together with evolutionary tracks from Ekstrom et al. 
(120121 ). Both figures also include the RSGs from our previous 
studies (VY CMa from Wittkowski et al. 120121 AH Sco, UY Set, 
KW Sgr from Arroyo-Torres et al. 2013 ). as well as Betelgeuse 
based on the data by Ohnaka et al. (120111) and VX Sgr based 
on the distance by Chen et al. (120071 ) and the angular radius by 
Chiavassa et al. (1201 Pal i, analyzed by us in Arroyo-Torres et al. 
(120131 ). Finally, we included the red giant stars from Marti-Vidal 
et al. (1201 lt > and Arroyo-Torres et al. (120141 ). All sources are con¬ 
sistent within their errors with the different calibrations of the 
effective temperature scale and with the red limits of the evo¬ 
lutionary tracks. The positions on the HR diagram of our new 
sources are close to evolutionary tracks corresponding to an ini¬ 
tial mass of 20-25 M 0 without rotation or 15-20 M 0 with rotation 
(V602 Car), 12-15M 0 without rotation or 9-15M e with rotation 
(HD 95687), 5-12M 0 without rotation or 7-9 M 0 with rotation 
(HD 183589). HD 183589 may thus also be a (super-)AGB star 
and not a RSG star (Siess l20l0l ). 

We note that the red giants with luminosities below 
log(L/L 0 ) ~ 4 are located systematically to the right of the 
Ekstrom tracks, and that a better agreement for these sources 
can be found using the STAREVOL grid (Lagarde et al. 
[2012) . which includes thermaline mixing unlike Ekstrom’s grid 
(Arroyo-Torres et al. 12014 ). 
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Table 4. Fundamental parameters of V602 Car, HD 95687, and HD 183589. 


Parameter 

V602 Car 

HD 95687 

HD 183589 

Ref. 

F bo , (10 1U Wm“-) 

11.30± 1.69 

4.84±0.73 

5.49±0.82 

1 

d (pc) 

1977±75 

1977±75 

621±62 

2 

L (10 31 W) 

5.28±0.89 

2.26±0.38 

0.25±0.06 

1,2 

log(L/L 0 ) 

5.14±0.17 

4.77±0.17 

3.82±0.25 

- 

Gross (mas) 

5.08±0.75 

3.17±0.50 

2.95±0.50 

This work 

R(Ro) 

1050±165 

674±109 

197±39 

2,4 

Teff (K) 

3432±263 

3467±282 

3709±322 

3,5 

iog(r eff ) 

3.54±0.08 

3.54±0.08 

3.57±0.09 

- 

log(g) 

-0.30±0.16 

-0.14±0.14 

0.80±0.17 

this work 

M(M 0 ) 

20-25 

12-15 

7-12 
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Notes. 1: Kharchenko ( f200TI . Cutri et al. 12003b . IRAS ( fl988l> ; 2: Humphreys et al. fl978i > - HD 95687, V602 Car (cluster CAR OB2); van 
Leeuwen i2007t - HD 183589; 6: Values obtain by the position of the stars in the HR diagram with the evolutionary tracks from Ekstrom et al. 
( 2012 ): We assumed a 15% error in the flux. The distance error was based on the values from Humphreys et al. ( 19781 ) by V602 Car and HD 95687. 
From HD 183589, we assumed a 10% error in the distance. The errors in the luminosity, effective temperature and radius were estimated by error 
propagation. 
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Fig. 10. Ratio between the square visibility in the continuum (av¬ 
erage between 2.27 and 2.28 pm) and the square visibility in 
the CO (2-0) line vs. log(L/L 0 ) (top) and the variability am¬ 
plitude (bottom) for a sample of RSGs in blue (Arroyo-Torres 
et al. 120131 Wittkowski et al. 120121 this work), AGB stars in 
red (Arroyo-Torres et al. 2014 ) and Mira variable stars in ma¬ 
genta (Wittkowski et al. 1 20111 ) for visibilities in the continuum 
between 0.2 and 0.4. The dotted black lines show the range of 
predictions by the best-fit PHOENIX model of these sources, the 
dashed orange the range of predictions by the best-fit RHD sim¬ 
ulation. Both lines overlap. 


4. Characterization of the extension of the 
molecular atmosphere 

We characterized the observed extensions of the molecular lay¬ 
ers in order to better understand how the fundamental parame¬ 
ters affect them. We also want to compare the behavior of RSG 
stars and Mira stars, which also exhibit extended molecular lay¬ 
ers (cf., e.g., Wittkowski et al. 1201 1 ). We used the ratio of the ob¬ 
served visibilities in the continuum band (average between 2.27 
and 2.28 pm) and the first (2-0) CO line at 2.294pm as an ob¬ 
servational indication of the contribution from extended atmo¬ 
spheric CO layers. Since this ratio depends on the value of the 
visibility in the continuum (V;r ont ), i.e. on how well the source 
was resolved, we limited the study to continuum squared visibil¬ 
ities between 0.2 and 0.4, a range where the visibility function 
is nearly linear. Although this approach may be limited by the 
limited spectral resolution of our observations of R ~ 1500 and 
by the low number of observations per source, it is appropriate 
for a first comparison of the extensions of RSG stars and other 
evolved stars. 

Fig.fTOlshows the resulting ratios (V" ont /V^ 0 ) for our sources 
vs. log(L/L 0 ) and AV, considering the RSGs by Wittkowski et 
al. (120121) . Arroyo-Torres et al. (120131) and this work (not repre¬ 
sent HD 183589 because their visibilities in the continuum are 
greater than 0.4), as well as the giants from Arroyo-Torres et 
al. (120141) . Results for Mira stars obtained by Wittkowski et al. 
(1201II) are also shown for comparison. Fig. [TT] shows the same 
visibility ratio but only for the RSG sample. In this case, we 
show the visibility ratio vs. L/L e , log(g), and 77 ir- Also shown 
are the ranges of the predicted visibility ratios based on the 
PHOENIX model atmospheres, as well as based on 3-D RHD sim¬ 
ulations, which will be discussed in the next section. 

Fig.flOlshows that the giant stars (red) are consistent with the 
PHOENIX models and the convection models (RHD simulations). 
[i Peg shows an atmospheric extension larger than predicted by 
the PHOENIX models, however it is relatively small and it does 
not show up with any significance using this metric of plotting 
(Vjoni/V ( 2 :0 ). The other giant stars do not show an extended CO 
band (the visibilities in the continuum are very similar to those 
in the CO (2-0) line). The RSG (blue) and Mira stars (magenta) 
show a more extended CO layer and therefore they are not con¬ 
sistent with hydrostatic models or convection models. The RSG 
and Mira stars show similar atmospheric extensions. 
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Fig. 11. Ratio between the square visibility in the continuum (av¬ 
erage between 2.27 and 2.28 yum) and the square visibility in the 
CO (2-0) line vs. log(L/L 0 ) (top), log(g) (middle), and the ef¬ 
fective temperature (bottom) for the sample of RSGs (Arroyo- 
Torres et al. 12013; Wittkowski et al. 120121 this work). The dotted 
black lines show the range of predictions by the best-fit PHOENIX 
model of these sources, and the dashed orange the range of pre¬ 
dictions by the best-fit RHD simulation. Both lines are overlap. 
The dashed gray lines, in the top and middle, show linear fits. 


Fig. HU suggests a correlation between the visibility ratio of 
our RSGs and the luminosity and surface gravity, which is indi¬ 
cated by linear fits (dashed gray lines). The correlations indicate 
an increasing atmospheric extension with increasing luminosity 
and with decreasing surface gravity. We do not observe a cor¬ 
relation between observed atmospheric extension and effective 
temperature or variability amplitude. In the case of Mira stars, 
we do not observe any correlation. The lack of a correlation be¬ 
tween visibility ratio and the luminosity and/or surface gravity 
suggests that different processes may be responsible for extend¬ 
ing the atmosphere in RSGs and Miras. On the other hand, we 
note that considerable atmospheric extensions for RSG stars are 


observed only for luminosities beyond ~ 1 * 10 5 L Q and for sur¬ 
face gravities below log(g)~0. 


5. Comparisons with convection and pulsation 
models 

In the following we discuss physical mechanisms that have been 
discussed as possible drivers of the observed large extensions of 
RSG atmospheres. In particular, we discuss convection models 
and pulsation models. Pulsation models have been successful to 
explain observed atmospheric extensions of Mira-variable AGB 
stars with shock fronts passing through the atmospheres. 

5.1. 3-D RHD simulations 

The 3-D radiation-hydrodynamics simulations of red supergiant 
stars were computed with C0 5 B0LD (Freytag et al. 12012) . The 
code solves the coupled equations of compressible (magneto- 
)hydrodynamics and non-local radiative energy transport on a 
Cartesian grid. The ”star-in-a-box” geometry is used and the 
computational domain is a cubic grid equidistant in all direc¬ 
tions; the same open boundary condition is employed for all 
sides of the computational box. The tabulated equation of state 
takes the ionization of hydrogen and helium and the formation 
of Hi molecules into account. The opacity tables are gray or use 
a frequency-binning scheme. 

The main model parameters are (Chiavassa et al. l201 lab : the 
stellar mass (entering the gravitational potential), the input lumi¬ 
nosity in the core, and the abundances that were used to construct 
the equation-of-state and opacity tables. Average values of stel¬ 
lar radius, effective temperature, and surface gravity have to be 
derived from a relaxed model (Chiavassa et al. 12009 , 201 la ). 

The models show large-scale convection cells (several in 
RSGs and only one to two for AGB stars) that span the entire 
convective envelope and have lifetimes of years. On the surface, 
there are smaller shorter-lived cells that resemble solar granules. 
According for instance to Samadi et al. (120011) . convective flows 
excite acoustic waves through turbulent Reynolds stresses (i.e., 
essentially velocities) and turbulent entropy fluctuations. In the 
part of the 3-D RHD simulations that comprises the stellar inte¬ 
rior, both are largest in or close to the downdrafts. Accordingly, 
acoustic waves are produced mostly in the downdrafts, and par¬ 
ticularly during merging of downdrafts. While the fast down- 
drafts themselves actually impede the outward propagation of 
the waves, the waves spread into the surrounding upflow regions 
in which they can easily reach the stellar surface. When these 
waves reach the thin and cold photospheric layers, they steepen 
into shocks that travel through the atmosphere with a decreasing 
density of the post-shock gas. Their induced dynamical pressure 
exceeds the thermal gas pressure leading to a significant increase 
in the density scale height. 

We used the pure-LTE radiative transfer OptimhD 
(Chiavassa et al. 12009 ) to compute intensity maps from 
three of these RHD simulations. This code take the Doppler 
shifts occurring due to convective motions into account. The 
radiative transfer equation is solved monochromatically using 
pre-tabulated extinction coefficients as a function of tempera¬ 
ture, density, and wavelength. The lookup tables are computed 
for the same chemical compositions as the RHD simulations 
using the same extensive atomic and molecular continuum and 
line opacity data as the latest generation of MARCS models 
(Gustafsson et al. 12008 ). 
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The three red supergiant (RSG) simulations are shown in the 
Table 0 On the other hand, we also tried the higher resolution 
simulation (401 3 grid points) st36gm00n05 (also in table 0. We 
considered several snapshots for each simulation and computed 
intensity maps in the wavelength range 1.90-2.60/rm with a con¬ 
stant spectral resolution of R-A/5 ,i ~ 20000. Moreover, for ev¬ 
ery wavelength, a top-hat filter including 5 wavelengths close by 
has been considered. In total, about 70000 images for every sim¬ 
ulation snapshot have been computed to cover the wavelength 
range of the observations. 

Then, using the method explained in detail in Chiavassa et 
al. (i2009t> . we computed azimuthally averaged intensity profiles. 
These profiles were constructed using rings regularly spaced in 
/r (where /j = cos{6) with 6 the angle between the line of sight 
and the radial direction). 

Finally, we averaged the monochromatic intensity profiles 
to match the spectral channels of the individual observations 
(we reduced the spectral resolution to AMBER observations). 
Afterward, we estimated the synthetic visibility of each baseline 
using the Hankel transform in the same way as for the PHOENIX 
models described above. 

Fig. El shows the intensity images at a continuum wave¬ 
length (2.20yum) and at the CO (2-0) line (2.294/rm), together 
with the azimuthally averaged intensity profiles for the example 
of the snapshot of model st35gm03n07. The intensity images il¬ 
lustrate that the intensity in the CO line is lower by a factor of 
about 2 compared to the intensity in the continuum, which is 
consistent with observed flux spectra such as in Fig. [T] for the 
example of V602 Car. In the images, the linear intensity range 
is between 0 and 130000 erg/s/cm 2 /A. Keeping the same scale 
seems to reduce the apparent luminosity of the CO image and 
thus also the contrast between bright and dark regions, however 
it is the opposite. The image of CO has a pseudo-continuum con¬ 
tribution due to CO opacities that increase the surface contrast 
and decrease the level of the structures, in particular at the limb. 
The CO line surface looks slightly more extended (purple color 
close to the stellar limb in the figure El central panel) but only 
by a few percent (~7%, estimated from Fig.[T2lbottom panel at 
the 0% intensity) 

Fig. El ( t0 P and middle panels) shows the example of our 
V602 Car data of Fig.|]]compared with the 3-D RHD simulation 
of Fig. El In the bottom panel, we show the predicted squared 
visibility amplitudes of the best-fit PHOENIX model to our data 
of V602 Car (cf. Fig. Q]i compared to those of the RHD simula¬ 
tion of Fig. El The model-predicted visibility curves of the 3- 
D RHD simulation are very similar to the hydrostatic PHOENIX 
model at the AMBER resolution and can thus not explain the 
large observed atmospheric extensions of RSG stars. 

The results are similar for all the RHD simulations employed 
in this work. In summary, it appears from the visibility curves 
that the atmospheric extension of RHD simulations is compa¬ 
rable to typical hydrostatic PHOENIX model within the spectral 
resolution of the actual observations. 

RHD simulations of RSGs do not lead to a considerable in¬ 
crease of the radius. However, the situation may be different for 
the AGB simulations (Freytag & Hofner [2008i ). where the ratio 
of the pressure scale height to the stellar radius is larger, giv¬ 
ing rise to relatively larger convective structures, to larger-scale 
waves and shocks, and finally to a noticeable increase of the stel¬ 
lar radius. Thus, shock waves in AGB simulations may explain 
the atmospheric extension we see in the observations of Mira 
stars. 


5.2. Pulsation models 

Self-excited pulsation models of Mira-variable AGB stars 
(Ireland et al. 120041 I2004bl 120081120TT1 Scholz et al. I20L41) 
have been successful to describe interferometric observations of 
these sources, including their extended atmospheric molecular 
layers (Woodruff et al. 120091 Wittkowski et al. 12011 . Hillen et 
al. 120121 . The observed visibility spectra of our RSG sources 
show similar features as those of Mira stars, in particular at the 
CO bandheads (see Sect. 3.3). We have thus as a first step inves¬ 
tigated whether pulsating model atmospheres of Mira variables 
can provide a good fit to our RSG stars as well, although the stel¬ 
lar parameters, in particular mass and luminosity, are very differ¬ 
ent and variability amplitudes of RSGs are typically lower by a 
factor of 2-3 compared to Miras (Wood et al. 1983 ). We used the 
recent CODEX model series by Ireland et al. ( 120081 1 2011 ) and, as 
an example, have found best-fit models to the data of V602 Car 
as shown in Fig. |T] Fig.fl4lshows one of the best-fit CODEX mod¬ 
els compared to these data. The shown model is model 261460 
of the o54 series. The o54 series is designed to describe the pro¬ 
totype Mira variable omi Ceti with a non-pulsating parent star of 
M=l.l M 0 , L=5400L o , R=216R e , P=330days. Model 261460 
is a model at phase 0.2 within a particularly extended cycle, with 
L-7420 L e and 7' c o -3154 K. The model does not show all ob¬ 
servational details, in particular the detailed overall slope of the 
flux and visibility spectra. However, it does match the drops of 
the CO bandheads in the visibility spectrum. Fig. fl5l shows the 
intensity profile of this model, in the continuum (2.25 /um) and in 
the CO (2-0) line (2.294 /um). as an illustration of which kind of 
atmospheric extension is required to describe the observed visi¬ 
bility data. The intensity profile in the CO line extends to beyond 
2 Rosseland radii, compared to an extension of a few percent as 
predicted by the 3-D RHD simulations of RSGs (Fig.EJ- 

As a next step, we calculated a new pulsation model based 
on stellar parameters that are typical for an RSG star: the non¬ 
pulsating parent star has M-15 M B , L=1.26T0 5 L Q , R=954R 0 , 
and its effective temperature is about 3600 K. The pulsation 
period of this model is about 750 days. These parameters are 
close to the star V602 Car discussed above (7)^=3432+263 K, 
L=( 1.35 +0.23)-10 5 L o , 7?=1050±1657? o , P=645days). Fig. El 
shows the radius variation of selected mass zones of this pul¬ 
sation model, the velocity at the photosphere, and the visual 
light curve. The amplitude of the photospheric radius variation 
is about 10% with radial velocities of up to about 5 km/sec. 
The model reproduces the amplitude of the visual light curve 
of V602 Car of about 1-1.5 mag. Whilst shock fronts enter the 
stellar atmosphere in a typical CODEX model of a Mira variable 
at or below optical depth l, leading to a geometric extension of 
the stellar atmosphere of the order a few Rosseland radii (e.g., 
Ireland et al. 1201T ). it turns out that no shock fronts reach at 
any phase the atmospheric layers in case of the RSG model. We 
also computed a higher amplitude version of this model with V 
mag., velocity, and radius amplitudes that are larger by a fac¬ 
tor of 2-3. Although the surface moves in and out more, there 
are still no shocks in the atmospheric layers. We note that the 
moving photosphere might induce shocks at higher radii, beyond 
the region of these models, where the density is very low. This, 
however, could also not explain the observed atmospheric exten¬ 
sions above the photospheric layers. In our test model, the ratio 
(Ktross = 10 4 ) - Ktross = 1))/Ktross = 1) which roughly 
measures the atmospheric extension is only 0.04 at four sam¬ 
ple phases (-0.02, 0.27, 0.52, 0.76), and we see hardly depth- 
dependent outflow/infall velocities below 4 km/s (Fig. El l. For 
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Table 5. RHD simulations of red supergiant stars used in this work. 


model 

grid 

[grid points] 

^pot 

[Mo] 

L 

[L 0 ] 

Teff 

[K] 

R* 

[R e ] 

logg 

[cgs] 

st35gm03n07" 

23S 3 

12 

91 932±1400 

3487±12 

830.0±2.0 

-0.335±0.002 

st35gm03nl3* 

235 3 

12 

89 477±857 

3430±8 

846.0±1.1 

-0.354±0.001 

st36gm00n06 

255 3 

6 

23 994±269 

3660±11 

391.3+1.5 

0.010±0.004 

st36gm00n05* 

401 3 

6 

24 233±535 

3710±20 

376.7±0.5 

0.047±0.001 


Notes, a) used in Chiavassa et al. l2009ll2010bl[201 lbl b) Chiavassa et al. 1201 la] 


comparison, the PHOENIX model plotted in Fig.Qjhas an exten¬ 
sion of 0.04 as well. 

A few test computations with higher parent-star luminos¬ 
ity (3 10 s L g ) and higher mass (25 and 15 M Q ) yield similar 
results. This means that the hydrostatic approximation should 
be well suited for calculating the stratification of the RSG at¬ 
mosphere. We note in this context that the atmospheric scale 
height, H oc T e $ * R 2 /M, which would be an approximate mea¬ 
sure of the atmospheric extension in case of a static stellar atmo¬ 
sphere, is of the same order of magnitude for typical Miras and 
for typical RSGs. For, e.g., the o54 model of the CODEX series 
that has about o Cet parameters and the above RSG model, one 
finds //(o54)///(RSG) = 0.7 and very small atmospheric exten¬ 
sions //(o54)/R(o54)~0.023//r and H(RSG)/R(RSG)~0.008/ju. 
The mean molecular weight is p ~ 1.3 for a normal solar element 
mixture with atomic hydrogen (molecular hydrogen is formed 
only in the high layers of the Mira models). We also note that 
H/R is approximately the ratio (gas particle thermal kinetic en¬ 
ergy )/(gravitational potential at the stellar surface). In the RSG 
model, the typical gas particle has a speed of ~8 km/s in the at¬ 
mosphere whereas pulsation velocities are only ~3 km/s. Hence, 
the pulsation energy is only -14% of the particle thermal energy 
and, therefore, does not lead to significant atmospheric extension 
of the models. 

Overall, the pulsation models for typical parameters of RSG 
stars lead to compact atmospheres with extensions similar to 
those of the PHOENIX and RHD models discussed above, and 
can thus also not explain the observed extensions of the molec¬ 
ular layers. 

5.3. Discussion on alternative mechanisms 

We showed in Sects. [3] and [4] that our interferometric observa¬ 
tions of RSGs provide evidence of extended molecular layers for 
a sample of RSGs with luminosities between about 1 x 10 s L 0 
and 5 x 10 5 L 0 and effective temperatures between about 3350 K 
and 3750 K. These extensions are comparable to those of Mira 
variable AGB stars, which typically reach to a few photospheric 
radii. Furthermore, our observations of RSGs indicate correla¬ 
tions of increasing atmospheric extension with increasing lumi¬ 
nosity and decreasing surface gravity, where considerable exten¬ 
sions of the CO (2-0) line are observed for luminosities above 
~ 1 * 10 5 L 0 and for surface gravities below log(g)~0. 

Comparisons of our interferometric data to available 1-D hy¬ 
drostatic PHOENIX model atmospheres, new 1-d pulsation mod¬ 
els, and available 3-D RHD simulations show that all of these 
theoretical models result in very similar synthetic visibility val¬ 
ues for the parameters of our observations. They indicate a com¬ 
pact atmospheric structure for all considered models, and none 


of these can currently explain the observed major extensions of 
the atmospheres of RSGs. 

A large number of spectroscopic and interferometric stud¬ 
ies of RGs and RGSs show that classical models do not explain 
the extended atmospheres of these sources (Tsuii 12003] 12008l >. 
Studies of resolved 12 /mi spectra of red giant and supergiant 
stars show strong absorption lines of OH and HiO, which are 
even larger than expected from a classical photosphere without 
a MOLecular sphere around the stars (Ryde et al. 12002), 2003 1 
I2006all2006bb . 

Hereby, the new 1-D pulsation test models show velocity am¬ 
plitudes below about 10 km/sec. These velocity amplitudes are 
consistent with observed long-term average velocity curves of 
Betelgeuse of ~ 9km/sec (Gray 12008,, Fig. 6) and of a sample 
of RSGs of < lOkm/sec (Josselin & Plez 12007 ). These 1-D pul¬ 
sation models are also consistent with typical amplitudes of the 
long-period visual light curves of RSGs of 1-3 mag. 

In addition to these long-term average velocity curves, which 
we interpret as being caused by ordered pulsation motions and 
which are well explained by the new 1-D pulsation models. 
Gray (120081 1 and Josselin & Plez (120071) also report on higher 
velocity gradients and turbulent motions in the atmospheres 
on short time scales with velocities of up to 30km/sec. These 
were explained by granulation and giant convection cells ac¬ 
companied by short-lived oscillations. In addition, Ohnaka et 
al. ( 20111 12013! ) reported on interferometric observations of 
Betelgeuse and Antares at two epochs each, which confirmed 
time variable atmospheric velocities of up to 20-30km/sec. 
They suggested that these motions are related to the wind driving 
mechanism. They also concluded that the density of the extended 
outer atmospheres of Antares and Betelgeuse are significantly 
higher than predicted by current 3-D RHD simulations, so that 
convection alone can not explain the formation of the extended 
atmospheres. This conclusion is confirmed by our direct com¬ 
parisons of interferometric data and synthetic visibilities based 
on current 3-D RHD simulations, and for a larger sample of 6 
additional RSGs. 

Current 3-D RHD simulations of RSGs are still limited in 
spatial resolution. However, we can compare older models of 
RSGs with low numerical resolution (235 3 grid points) and new 
models with better resolution (401 3 grid points). On the other 
hand, we also have better resolved models of AGB stars and 
even better resolved local models of solar-type stars - with a 
range of resolutions. From these comparisons, we conclude that 
future RSG simulations with higher resolution are indeed im¬ 
portant and desirable, but that we do not expect the atmospheric 
velocity fields due to convection and pulsations alone to grow 
enough to remotely reach the amplitude necessary to give a mol- 
sphere extension of a few stellar radii as observed. 
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Josselin & Plez d2007l > suggested that high velocities and 
steep velocity gradients, possibly caused by convective motion, 
generate line asymmetries, that turbulent pressure decreases the 
effective gravity, and that this decrease combined with radiative 
pressure on lines initiates the mass loss. The spherically symmet¬ 
ric CODEX Mira model atmospheres do, indeed, show noticeable 
radiative acceleration a rM \ which significantly affects the effec¬ 
tive gravity in the atmospheric layers (e.g.. Fig. 6 with Eq. 1 in 
Ireland et al. 12008 1. but still without any net outward acceler¬ 
ation. One should realize that the pulsation-generated very ir¬ 
regular velocity stratification (resulting in irregular p(r), e.g.. 
Figs. 14/15/16 in Ireland et al. [20111 and Fig. 2 in Scholz et 
al. 120141 ) favors large a rac j because of substantial line shifts. 
Velocity (outward/inward) variations can be quite strong, up to 
a few km/s between shock fronts, apart from the big “jumps” 
at shock front positions. The new RSG pulsation models de¬ 
scribed in Sect. E 2 ] do not show this velocity stratification be¬ 
havior. However, they are also solely based on long term velocity 
variations caused by pulsation motion. It is plausible that much 
steeper velocity gradients on shorter time scales, as reported by 
Gray (120081 ). and not included in these models, might generate 
accelerations on Doppler-shifted molecular lines that sufficiently 
extend the molecular atmosphere. Although already suggested 
in 2007, more detailed calculations of this effect for RSGs and 
its implementation in the 3-D RHD simulations are still pending 
and are also outside the scope of the present paper. Our observed 
correlation of increasing atmospheric extension with increasing 
luminosity and decreasing surface gravity (cf. Fig. ITU supports 
such a radiatively driven extension of the atmospheres of RSGs. 

The proposed scenario of radiation pressure on Doppler- 
shifted lines is in a a way reminiscent of what happens in the 
winds of hot stars. However, it should be noted that -unlike for 
hot stars- RSGs form dust at larger radii (typically ~20 stellar 
radii, Danchi et al. 1994 1 so radiation pressure on dust cannot 
occur in the wind acceleration zone (Josselin & Plez, 1 2007 ). In 
the case of RSGs, the proposed radiative pressure on molec¬ 
ular lines may only levitate the molecular atmosphere up to 
radii where dust can form (as suggested by Josselin & Plez, 
12007 ). analogous to shock fronts for Miras, albeit not lifting 
the material outside the gravitational potential. However, the 
detailed dust formation, condensation sequence, and mass-loss 
mechanisms are also not yet fully understood in the case of 
oxygen-rich AGB stars (cf., e.g., Karovicova et al. 120131 for a 
recent discussion). Recent indications based on a polarimetric 
aperture masking technique (Norris et al. 120121 ) and on mid- 
infrared interferometry (Wittkowski et al. 2007 ; Karovicova et 
al. 1201 ll 120131 ) point toward transparent dust grains forming al¬ 
ready at relatively small radii of about 1.5 stellar radii. These 
may be grains of amorphous ALO 3 and/or magnesium-rich 
iron-free (“forsterite”) silicates, while iron-rich silicates form 
at larger radii. Verhoelst et al. (120061) found a similar evidence 
for amorphous alumina in the extended atmosphere of the RSG 
Betelgeuse as close as ~ 1.4 stellar radii, while iron-rich silicates 
indeed had an inner radius of 20 stellar radii in their model. 
However, Kaminski et al. (120131 ) do not observe the presence 
of AlO in the inner outflow in the spectrum obtained with the 
Ultraviolet and Visual Echelle Spectrograph (UVE) at the Very 
Large Telescope (VLT). On the other hand, they have detected 
AlO-bearing gas in the wind-acceleration zone, out to 20 R„. If 
the presence of transparent grains is confirmed at small radii of 
~ 1.5 stellar radii, it would play a crucial role for the dust con¬ 
densation sequence and the overall mass-loss process of both 
oxygen-rich AGB stars and RSGs. 


Grunhut et al. (120101) . Auriere et al. (1201 Ok and Bedecarrax 
et al. (120131 ) detected weak (<1G) magnetic fields in one third of 
their sample of late type supergiants, with their sensitivity limit. 
They suggested that magnetic fields are generated by dynamo 
action, and that weak fields may be excited in all cool super¬ 
giants. It has been discussed that magnetic fields may also play 
a role in the mass-loss process of RSG stars alongside other pro¬ 
cesses as discussed above. Thirumalai & Heyl (120121) recently 
presented a magnetized hybrid wind model for Betelgeuse, com¬ 
bining a Weber-Davis stellar wind with dust grains, that was able 
to lift stellar material up from the photosphere and into the cir- 
cumstellar envelope. One of us (BF) recently produced a set of 
5 -Mq-RSG RHD models without and with magnetic field. These 
first simulations with magnetic fields did not significantly alter 
the stratification of the model, however, more calculations are 
need to reach a real conclusion. We also note that the good over¬ 
all fit of the hydrostatic PHOENIX models in the continuum 
bands means that the temperature structure vs. optical depth is 
roughly correct, so that additional heating of the atmosphere is 
not a missing mechanism. In addition, rigid rotation or differen¬ 
tial rotation (e.g., a rapidly rotating core) are processes that are 
not included in our models but may contribute to the atmospheric 
extension and to the mass-loss process. 

6. Summary and conclusions 

Our spectro-interferometric observations with AMBER are a 
good tool for studying the continuum-forming layer and molec¬ 
ular layers separately and for constraining the atmospheric struc¬ 
ture of RSGs. We used the continuum near 2.2pm, which is 
mostly free of molecular contamination, to estimate the angular 
diameters of the targets. Together with estimates of the distances 
and the bolometric fluxes, we also derived fundamental stellar 
parameters such as the luminosity, the Rosseland radius, and the 
effective temperature. 

With the effective temperature and the luminosity, we lo¬ 
cated our targets in the HR diagram. Their locations are close 
to evolutionary tracks that correspond to initial masses of 20 - 
25/15-2OM 0 (V6O2 Car), 12-15/9-15M 0 (HD 95687), and 5- 
12/7-9 M 0 (HD 183589) with or without rotation. All target po¬ 
sitions are consistent with the red limits of recent evolutionary 
tracks. HD 183589 shows a lower luminosity and thus lower 
mass compared to our other sources. It may more likely be a 
(super-)AGB star instead of a RSG star. 

The near-infrared spectra of all our target are reproduced 
well by hydrostatic PHOENIX model atmospheres, including the 
CO bands. We had observed the same behavior in our previous 
project (Wittkowski et al. 12012 1. Arroyo-Torres et al.l 2013 ). 

The observed visibility curves of our sample of RSGs show 
large drops in the CO (2.3-2.5 pm) and partly in the H 2 O (1.9- 
2.1pm) bands, indicating major extensions of the atmospheric 
molecular layers. As a first characterization of the extensions, we 
calculated the square visibility ratios between the nearby contin¬ 
uum and the CO (2-0) bandhead. More detailed observations us¬ 
ing a larger number of data points, possibly imaging campaigns, 
and a higher spectral resolution to better isolate the CO band- 
head are desirable for a more accurate description of the exten¬ 
sions. Nevertheless, we observed a linear correlation between 
the visibility ratios of our RSGs and the luminosity and surface 
gravity, indicating an increasing atmospheric extension with in¬ 
creasing luminosity and decreasing surface gravity. With that, 
we observed considerable atmospheric extensions of RSGs only 
for luminosities beyond ~ 10 s L 0 and for surface gravities below 
log(g) ~0. We did not observe a correlation with effective tem- 
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perature or variability amplitude. Compared to Mira stars, we 
noticed comparable extensions between RSGs and Mira stars, 
which extend to a few stellar radii. Mira stars did not exhibit 
the correlations with luminosity and surface gravity that we ob¬ 
served for RSGs. This suggests that the extended atmospheric 
structure may be generated by different processes for each type 
of star. For Miras, the extensions are believed to be triggered by 
shocks that are generated by pulsations and that enter the atmo¬ 
spheric layers. 

The synthetic visibility amplitudes of hydrostatic PHOENIX 
models did not predict the strong visibility drops in the molecu¬ 
lar bands, and can thus not explain their observed major exten¬ 
sions. To further constrain processes that were discussed as pos¬ 
sible drivers of the extensions, we compared our data to avail¬ 
able 3-D RHD simulations of convection and to new 1-D pulsa¬ 
tion models for typical parameters of our RSGs. Both models re¬ 
sulted in a compact atmospheric structure as well, produced sim¬ 
ilar observable synthetic visibility values as the PHOENIX mod¬ 
els, and could therefore not explain the major extensions. We 
note that the actual atmospheric extension of 3-D RHD simula¬ 
tions is not enough to take the contribution of the molecular ex¬ 
tended layers into account. Improvements in this direction may 
alter the appearance of the stellar surface with respect to typical 
3-D RHD simulations of the photosphere as shown in Fig. [T2] 
This may be particularly important for interferometric observa¬ 
tions taken with broad bandpasses, which are contaminated by 
molecular bands. 

Our observed correlation of increasing atmospheric exten¬ 
sion with increasing luminosity and decreasing surface gravity 
supports a scenario of a radiatively driven extension caused by 
radiation pressure on Doppler-shifted molecular lines. We spec¬ 
ulate that another ingredient of the mass-loss process could be 
acceleration on dust grains that may form already at a few stel¬ 
lar radii, a process which may possibly be further supported by 
magnetic fields or, differential rotation (e.g., a rapidly rotating 
core). In this case, the radiative acceleration on molecular lines 
would only be needed to levitate the atmosphere up to the point 
where the dust grains are formed. 
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Fig. 12. 3-D radiative-hydrodynamic al (RHD) simulations of 
convection for one snapshot of model st35gm03n07 (see table 
|5}. (Top) Image of the intensity at a continuum wavelength of 
2.20gm. (Middle) Image of the intensity at the CO (2-0) line at 
2.294/an. (Bottom) Azimuthally averaged intensity profiles of 
both images, where the black line denotes the continuum wave¬ 
length and the red line denotes the CO line. The radii is defined 
by r/R t = ^Jl - /u 2 (// is explained in the text). 
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Fig. 13. Observed normalized flux (top), square visibility ampli¬ 
tudes (middle) of V602 Car as in Fig. |T] compared to the predic¬ 
tion by the 3-D RHD simulations of Fig. [12] The black curves 
are the observational data, the green curve the RHD simulation, 
and the blue curve the best-fit UD model (it is underplotted to 
the green curve). Bottom: Enlargement of the synthetic squared 
visibility amplitudes of the CO (2-0) line (2.294/an) predicted 
by the PHOENIX model atmosphere that best fits our V602 Car 
data (cf. Fig. [T]) compared to the prediction by the RHD simula¬ 
tion of Fig. [T2] The stellar parameters of the PHOENIX mode are 
r eff =3400K, log(g) =-0.5, M=20M q . Those of the RHD simu¬ 
lation are r e ff=3487 K,; log(g)=-0.335, R=830R G , M —12 M 0 ). 
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Fig. 14. Observed normalized flux (top) and squared visibility 
amplitudes (bottom) of our V602 Car from Fig. [T] compared to 
one of the best-fit CODEX model atmospheres. The black lines 
denote the observational data, and the blue lines the model pre¬ 
dictions for all three baselines (bottom). The model parameters 
are listed in the main text. The thick vertical lines indicate the 
uncertainties at three wavelength intervals. 
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Fig. 16. Pulsation model of a RSG with M=15 M 0 , 
L= 126000 L®, Teff~3600K. Bottom panel: Radius variation 
of selected mass zones in a pulsating supergiant model with 
M = 15 M e and L = 126000 L 0 (black curves). The red curve 
is the position of the photosphere (defined as the layer where 
the Rosseland optical depth equals 2/3). Middle panel: The ve¬ 
locity at the photosphere. Top panel: The visual light curve of 
the model, where the bolometric correction is obtained from the 
tables in Houdashelt et al. ( 20QOal[2000bl . 


Fig. 15. Intensity profile of the CODEX model o54/261460 in the 
continuum at 2.25 fim (solid blue line) and at the CO (2-0) line 
at 2.294yum (dashed red line). 
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Fig. 2. As Fig. [U but for data of V602 Car obtained with the MR-K 2.1 pm setting on 4 April 2013 (left) and with the MR-K 2.3 pm 
setting obtained on 4 April 2013 (right). 
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Fig. 3. As Fig. |TJ but for data of HD 95687 obtained with the MR-K 2.1 pm setting on 4 April 2013 (left) and with the MR-K 2.3 
pm setting on 4 April 2013 (right). 
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Fig. 4. As Fig.Q] but for data of HD 95687 obtained with the MR-K2.3 pm setting on 25 April 2013. 
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Fig.5. As Fig.Q] but for data of HD 183889 obtained with the MR-K 2.1 pm setting on 4 August 2013. 
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Fig. 6. As Fig.E but for data of HD 183589 obtained with the 
pm setting on 29 July 2013 (right). 
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